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We present an extension of the velocity-dependent one-scale model for cosmic string evolution, 
which is suitable for describing the evolution of local and global monopole networks. We discuss the 
key dynamical features that need to be accounted for, in particular the fact that the driving force 
is due to the other monopoles (rather than being due to local curvature as in the case of extended 
objects) and new forms of energy loss terms due to monopole-antimonopole capture and annihilation. 
For the case of local monopoles we recover and generalize the results of Preskill, suggesting that the 
scaling law for the monopole correlation length is very sensitive to the annihilation rate. On the other 
hand, for global monopoles the long-range forces generically lead to linear scaling (just like in the 
case of local cosmic strings) . In this case we also find good qualitative agreement between our results 
and the numerical simulations of Bennett & Rhie and Yamaguchi, although future high-resolution 
simulations will be needed for quantitative comparisons. 



I. INTRODUCTION 

Since the pioneering work of Kibble [l| it has been 
known that topological defects necessarily form at cosmo- 
logical phase transitions if they are stable. The details of 
the phase transition, in particular the specific symmetry 
being broken, will determine what kind of defect net- 
work is formed, and whether or not the network is long- 
lived. Understanding defect formation, dynamical prop- 
erties and evolution is therefore a key part of any effort 
to understand the early universe. A thorough overview 
of the subject can be found in the book by Vilenkin and 
Shellard Q. 

In the past three decades, most of the work on defects 
has focused on cosmic strings, for the good reason that 
they are usually cosmologically benign and are a generic 
prediction of inflationary models based on Grand Unified 
Theories [!, 0| or branes [H, 0]. On the other hand, do- 
main walls and monopoles are cosmologically quite dan- 
gerous, so any models in which they arise are subject 
to very tight bounds. Nevertheless, understanding their 
detailed properties is still important. On the one hand, 
in order to confidently impose bounds on models which 
form them one must know in detail how they evolve. On 
the other hand it is becoming increasingly clear, partic- 
ularly in the context of models with extra dimensions 
such as brane inflation, that hybrid defect networks will 
often be produced. Two examples are semilocal strings 
and cosmic necklaces 0, H, H|- To understand the evo- 
lution of these hybrid systems one needs to understand 
not only the dynamics of the strings but also that of the 
monopoles. This article is a first contribution to this 
task. 
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The existence of magnetic monopoles was first sug- 
gested by Dirac in 1931 [lOj], but first demonstrated in a 
field theory context only in 1974 by 't Hooft and Polyakov 
[Til [l2| . It soon became apparent that they are cosmo- 
logically disastrous, and indeed this has beengiven the 
status of a 'problem' — the monopole problem |2l.ll3l| . For 
this reason, their study has been comparatively neglected 
relative to that of cosmic strings and even domain walls. 

Global monopoles have logarithmically divergent en- 
ergy due to the slow fall-off of angular gradients, which 
makes some of their properties counterintuitive. Even 
their stability was the subject of some debate because a 
global monopole with its core position fixed has an an- 
gular zero mode [l4[ and a finite energy barrier between 
topologically distinct sectors In reality, any angu- 

lar deformation away from spherical symmetry leads to a 
displacement of the core and the only decay channel left 
is monopole annihilation [l6|, [l7 1. Global monopoles are 
known to have scaling solutions H, [lj| , and could play a 
role in structure formation that can be constrained with 
CMB data [M [H| (see also H3). 

In this paper we propose an extension of the velocity- 
dependent one-scale model for string evolution [H, [24[ 
that can be used to study the evolution of local or global 
monopole networks. In future work we shall discuss how 
this can be applied to the study of hybrid defect net- 
works. 

The rest of the article is organized as follows. We start 
in Sect. |TT] by presenting a simplified derivation of the 
evolution equations for the two dynamical quantities of 
our model, a length scale that can be thought of as a cor- 
relation length and a root-mean squared (RMS) velocity. 
We then focus on the specific case of monopoles, first 
by providing (in Sect. IIIip a brief overview of previous 
work on monopole properties and dynamics and then by 
discussing (in Sect. IIV[) how these features can be en- 
capsulated in our model. In Sect. |V]we use the model 
to discuss the evolution of networks of local monopoles; 
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we shall see that our results agree with but generalize 
those of Preskill [25|. Similarly, the evolution of global 
monopole networks is discussed in Sect. IVII here we find 
good qualitative agreement with previous numerical sim- 
ulations by Bennett and Rhie [l8| and by Yamaguchi 
[l9l ] , although the relatively low resolution of these sim- 
ulations does not allow quantitative comparisons (which 
must be left for future work). In Sect. IVIII we revisit 
some our modeling assumptions on the energy loss terms 
and driving forces, and discuss how these affect our re- 
sults. Finally in Sect. I Villi we summarize our results and 
discuss some open issues. 



II. ANALYTIC MODELING GENERALITIES 

Broadly speaking, the key idea behind the analytic 
modeling of defect networks is that one give up the idea 
of studying their 'statistical physics' (which to a large 
extent can only be done numerically) and instead con- 
centrates on its 'thermodynamics'. In other works, one 
or more macroscopic quantities are chosen to describe 
the network, and the knowledge of the microphysics is 
used to derive evolution equation for these macroscopic 
quantities. This has the advantage of leading to rela- 
tively simple models which can in principle encapsulate 
most of the relevant physics, but it also has an associated 
cost: in going from the microphysics to the macrophysics 
one is forced to introduce phenomenological parameters, 
which can only be calibrated by direct comparison with 
numerical simulations. 

The first effort along these lines is Kibble's one-scale 
model of string networks [26[ (extended in J27[), which 
has a single macroscopic parameter: a lengthscale that 
can be identified as the string correlation length, the 
string curvature radius, or the inter-string distance — in 
the model's approximation they should be seen as iden- 
tical or at least comparable. This was later generalized 
to a three-scale model by Austin, Copeland and Kib- 
ble [28[ where there are three distinct lengthscales: the 
first two are again the string correlation length and the 
inter-string distance, while the third is a typical scale for 
small-scale 'wiggles'. 

A different approach recognizes that in order to be 
able to quantitatively describe the whole cosmological 
history of these networks one must be able to describe 
the evolution of the defect velocities, not least because 
depending on the cosmological epoch and on their own 
properties, the defects may be moving at non-relativistic 
or at ultra-relativistic speeds. This is the basis for the 
velocity-dependent one-scale model of Martins and Shel- 
lard [23. l24|. which retains Kibble's assumptions on the 
existence of a single lengthscale but adds the RMS ve- 
locity as a second macroscopic quantity. This model can 
be accurately derived starting from the Goto-Nambu ac- 
tion, and it has also been successfully tested against high- 
resolution field theory [2^| and Goto-Nambu numerical 
simulations [13, HH . More recently, it has also been ex- 



tended to domain walls 32]. In what follows we provide 
a simplified derivation of the dynamical equations of the 
model, along the lines of [13] . 

Consider a network of defects with n-dimensional 
worldshccts (n — 1 for monopoles) evolving in (3 + 1) 
space-time dimensions. Temporarily assume them to 
have velocity v, to be non-interacting and (for extended 
objects) planar. Then the momentum per unit comoving 
defect volume -simply the momentum, for monopoles- 
goes as 



p cx a =>• v 7 oc a 
from which we get by differentiation 
dv 



dt 



+ nH(l - v 2 )v = 0. 



(1) 



(2) 



On the other hand, under the above hypotheses the aver- 
age number of defects in a fixed comoving volume should 
be conserved, which implies 



p oc ja 



-(4-n) 



(3) 



and again, differentiating and using the velocity equation, 
we get 



^ + H[(4 - n) + nv 2 ]p = . 



(4) 



The hypotheses so far are, of course, widely unrealistic. 
However, we can use this as a starting point to build a 
reasonable model. As has been pointed out above, the 
validity of this process can be checked for the case of 
cosmic strings, where a more rigorous derivation has been 
done. 

Let us start by defining a characteristic length scale 



(•5) 



where M will have dimensions appropriate for the defect 
in question (i.e., monopole mass, string mass per unit 
length, or wall mass per unit area), and can also be writ- 
ten 



(6) 



with r\ being the symmetry breaking scale. Also, we in- 
terpret the velocity as being the RMS velocity of the de- 
fect network, and allow for energy losses due to interac- 
tions, which for extended defects can usually be modeled 
(purely on dimensional grounds) by 



dp v 



(7) 



we shall show below that this p dependence also applies 
to global monopoles while for local monopoles the energy 
loss is proportional to p 2 . More importantly, defects will 
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be slowed down by friction due to particle scattering. 
This can be characterized by a friction length scale 



M 



where we are defining 



M 



cx a 



Tt + l 



(8) 



(9) 



and 9 is a parameter counting the number of particle 
species (or degrees of freedom) interacting with the de- 
fect. 

We can also define an overall damping length which 
includes both the effect of Hubble damping and that fric- 
tion due to particle scattering 



rr nH+ T, 



(10) 



It is important to compare the relative importance of the 
two effects. Since the friction length scale will in most 
circumstances grow faster than the Hubble length, it is 
expected that friction will be dominant at early times, 
while Hubble damping will dominate at sufficiently late 
times. This is easy to confirm. Using the relations [l3} 



t 



for the radiation era, and 



mpi 

T 2 



(11) 



if in the radiation era, otherwise at 

T / 77 \^ ( T eq 



mpi 



mpi 



mpi 



(16) 



in the matter era. 

Putting together all of the above effects, we find the 
following evolution equation for the characteristic length 
scale L and RMS velocity v 

(4-n)^ = (4-n)HL + v 2 ^- + cv (17) 
at id 



(18) 



where in the latter we have included the possibility of 
further driving forces affecting the defect dynamics. Note 
that / has the units of acceleration, ic it is the force per 
unit mass. For extended objects (walls and strings) that 
have been extensively studied in the past, this driving 
force is obviously the local curvature, and we have 



1 L 



(19) 



we are implicitly assuming that our characteristic length 
scale is the same as the defect curvature radius. For 
monopoles the situation is more complicated, since there 
are forces due to other monopoles. This will be discussed 
below. 



mpi 



(12) 



for the matter era, we find that the ratio of the Hubble 
and friction length scale is given by 



1 



/ mpi 

l f H ~ ^ V 



in the radiation era, and 



T 



T 

mpi 



ji— 1/2 
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1/2 



(13) 
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in the matter era, where again r\ is the scale of the phase 
transition producing the defects. In the interest of sim- 
plicity, these estimates neglect couplings as well as num- 
bers of interacting species (which may not necessarily 
be of order unity) . For defects formed at the GUT phase 
transition, the friction length scale is ~ 10 3 times smaller 
than the Hubble length. Note that the ratio is constant 
for monopoles in the radiation era - hence monopoles will 
always be friction-dominated in the radiation era (but not 
in the matter era). Apart from this special case, friction 
domination ends at 



T 

mpi 



V 
mpi 



(15) 



III. OVERVIEW OF MONOPOLES 

We now provide a brief overview the properties of local 
and global monopoles. This is by no means exhaustive — 
we shall only focus on the dynamical properties that will 
turn out to be relevant for our analytic model. We shall 
follow the more extensive reviews in 0, [H, 34] , and refer 
the interested reader to them for a more detailed discus- 
sion. 

The local monopole solution has two characteristic 
length scales, r s and r v , identifying the radii of regions 
in which the scalar and vector fields depart significantly 
from their asymptotic behaviour. The monopole mass is 
approximately given by 



4?r 



47T 



->ln 



(20) 



Their initial separation, can range from the monopole 
thickness (for a second-order phase transition) to the 
horizon size (for a strongly first-order phase transition). 
The expansion of the universe will dilute the monopoles, 
so naively we might expect a solution of the form £ cx 
T _1 cx a. 

The force between local monopoles is just electromag- 
netic attraction, 



Femg (?") 



9 



(21) 
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where the magnetic charge is given by g — An/e. A 
monopole moving through a plasma also experiences a 
drag force due to its interaction with charged particles; 
the corresponding force is 



plasma 



(22) 



in agreement with Eqn. ([8]). The initial monopole den- 
sity is tim oc £~ 3 , though note that due to the Coulomb 
forces the positions of monopoles and antimonopoles are 
strongly anti-correlated: except, possibly, for a short 
transient, the nearest neighbour to a monopole is likely 
to be an anti-monopole. The Kibble mechanism would 
naturally produce such an anticorrelation already at for- 
mation. On the other hand, for slow transitions there is 
a competing mechanism in which the monopoles are the 
result of thermal fluctuations of the gauge field [H, [36| 
and they form in same-charge clusters. However, even 
in this case the typical number of monopoles per clus- 
ter is of order unity at the relevant energy scales around 
the GUT scale [36[. Then the attractive forces between 
monopoles and antimonopoles will dynamically establish 
the anti-correlations very quickly. 

One might naively expect monopole motion in a 
plasma to be like Brownian motion of heavy dust par- 
ticles in a gas or liquid. If this is so, then they would 
typically move with thermal velocities 



T 

v T ~ i — 

m 



1/2 



with mean-free path 



— (-) 

6T\t) 



1/2 



(23) 



(24) 



However, the Coulomb forces will again introduce a bias 
in the random walks, and so the defects gradually drift 
towards each other. The average drift velocity may be 
simply estimated from the balance between the electro- 
magnetic and drag forces, leading to 



(25) 



Q T 2 r 2 



This dissipates some energy and leads to the formation 
of bound monopole-antimonopole states, at a capture ra- 
dius 



T 



(26) 



Once the bound pair is formed, the monopole and anti- 
monopole spiral in, losing energy to the plasma drag and 
radiation, and eventually annihilate into gauge bosons. 
The rate of radiative energy loss into gauge bosons can 
be estimated using the classical electromagnetism radia- 
tion formula 



M 2 

67T 



(27) 



(a being the monopole acceleration), which leads to 



The lifetime of the bound state is 



9" 



9 2 



(28) 



(29) 



Notice that radiation losses increase strongly as the sepa- 
ration decreases, and indeed almost all the energy is lost 
after capture. 

These bound states are generically short-lived, al- 
though if desired one could design models where the an- 
nihilation would occur much later than the capture. In 
any case, this subtlety can in practice be neglected, since 
a captured pair is doomed to annihilate anyway — it can 
therefore be considered as decoupled from the rest of the 
network from the moment of capture onwards. 

This, together with the fact that most radiation losses 
occur in bound pairs, means that we need not model 
these effects explicitly. All we require is an adequate ac- 
count of the 'pair losses', that is we need to keep track of 
the rate at which bound pairs are forming and leaving the 
'free' monopole network. This is analogous to neglecting 
the effect of the small string loops when modeling the dy- 
namics of a network of 'infinite' strings — there one also 
does not need to explicitly account for losses to gravita- 
tional radiation (which occur mostly in the loops rather 
than in the infinite strings). 

Finally, note that the diffusive capture is only effective 
if the monopole mean free path in the plasma is smaller 
than the capture radius, which corresponds to 



T > 



6>V 6 2 g 3 



(30) 



At lower temperatures the monopoles and antimonopoles 
can only capture each other by emitting radiation. For 
thermal incident velocities we have already seen that the 
energy losses are very small, so the annihilation time is 
much larger than the expansion rate and the monopole 
to entropy ratio becomes effectively frozen. 

More interesting still is the case of global monopoles. 
The total energy of a monopole out to a given distance 
grows linearly with distance, so in some sense global 
monopoles are similar to local strings. The energy of 
a monopole-antimonopole pair at a distance R is 



£ 



AirrfR , 



and the corresponding attractive force is therefore 

d£ 



global 



OR 



(31) 



(32) 



which is independent of distance. This is analogous to 
monopoles connected by strings. But unlike that case, a 
global monopole is not paired with any particular anti- 
monopole, and so it is not a priori clear how efficiently it 
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can find a partner. Still, the existence of long-range forces 
means that capture and annihilation will be more effi- 
cient than in the local case, and therefore no 'monopole 
problem' is expected here. The rate of energy loss of 
a monopole-antimonopole bound state into Goldstone 
bosons has been estimated to be 1371 



e, 



gold 



->r 



(33) 



and its lifetime is t ~ R. The pair is expected to move at 
ultra-relativistic speeds and annihilate within a Hubble 
time. Again, this allows us to model the energy losses as 
going into bound pairs as opposed to modeling radiation 
losses explicitly. 

Finally, let us note that there is some previous nu- 
merical work on global monopole networks. Numerical 
simulations in the non-linear sigma model approxima- 
tion by Bennett and Rhie [l8j suggest that a scaling solu- 
tion is reached with a few monopoles per horizon volume, 
namely 



n H d 3 H = 3.5 ± 1.5 



(34) 



in the radiation era and 



n H d 3 H = 4.0 ± 1.5 (35) 

in the matter era. More recently, Yamaguchi [l!| gives 
for the quantity £ = nt 3 , 



£ r = 0.43 ±0.07 



(36) 



£ m = 0.25 ± 0.05 (37) 
and also attempts to measure velocities, finding 

u r = 1.0 ±0.3 (38) 



v m = 0.8 ± 0.3 . 
Note that for a scale factor a(t) oc t x 

d H {t) 



t 



1-A 



(39) 



(40) 



so the two sets of simulations agree with each other re- 
markably well in the radiation era, and broadly so in the 
matter era. 



IV. MODELING MONOPOLES 

We are now in a position to bring together the results 
of the two previous sections and discuss how the velocity- 
dependent one-scale model can be extended to describe 
the evolution of monopole networks. From the above 
discussion, it is clear that the force between a pair of 
local monopoles has the following form 



f local 



rjL 2 



(41) 



For global monopoles the situation is slightly more subtle, 
but also more interesting. The force between a pair of 
them is independent of distance, but recall that their 
mass grows proportionally to the distance. Therefore the 
acceleration is in fact inversely proportional to distance, 



f global 



k 

L 



(42) 



Again we see that they are in some sense like local strings. 
This is not too difficult to understand: a monopole will 
be effectively heavier when seen on larger scales, and its 
acceleration should therefore be correspondingly smaller. 

The next issue is the fact that there can be many 
monopoles and anti-monopoles in a given Hubble volume, 
so the various forces acting on a given one will partially 
cancel each other. A simple way in which one can try to 
model this as a 1/V~N effect. In other words, the accel- 
eration / becomes / / y/N, where the number of defects 
N in a Hubble volume d 3 H is given by 



(43) 



This should be fine for global monopoles (although it is 
an issue that warrants testing in numerical simulations). 
On the other hand, in the local case the situation is again 
more subtle due to the existence of anti-correlations in 
the positions of monopoles and anti-monopoles. The 
analysis of [38| indicates that the number of defects in 
that case is instead approximately given by 



Ni 



dn 
L 



(44) 



This is to be expected: since the nearest neighbour to 
a monopole is likely to be an antimonopole (and vice- 
versa) , then typically the attractive forces between neigh- 
bouring pairs will be larger than in the uncorrelated case. 
In other words, the cancellation mechanism is less strong, 
which is equivalent to saying that the effective number 
of neighbours is smaller. 

Finally, there is the issue of energy losses due to 
monopole annihilations. The generic form given by Eqn. 
(O, a natural extension of those that can be derived 
for strings and domain walls, is valid for the case of 
global monopoles (again, these are in some sense like lo- 
cal strings). For a single monopole-antimonopole pair 
p/p ~ £/£, and Eqns. (|3"Tj) and (|3"3"|) lead to p cx —p/R. 
This is another way of saying that the timescale for en- 
ergy losses corresponds (in the fundamental units we are 
using) to the lenghtscale R. On the assumption that the 
two lengthscales are comparable R ~ L and this matches 
Eqn. ([7]) apart from the allowance for generic velocities. 
Note that here we are assuming that the separation be- 
tween monopole-antimonopole pairs is comparable to the 
network's characteristic lengthscale. Although the two 
need not be the same, this is a valid assumption in the 
context of the simple one-scale model we are considering, 
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and additionally any discrepancy could to some extent be 
absorbed by a redefinition of the numerical coefficient c. 

On the other hand, in the local case the Coulomb forces 
between the monopoles and antimonopoles lead to a dif- 
ferent en ergy loss rate. Early work of Zel'dovich and 
Khlopov [39( | assuming that the monopole-antimonopolc 
annihilations were purely determined by the Coulomb at- 
traction of the magnetic charges found that the number 
density Um should evolve as 



dn M 
dt 



A T 3 > 



(45) 



a subsequent more detailed study by Preskill [25| , which 
implicitly allows for the effect of the anti-correlations, 
leads to the more generic form 



dn 



M 



<ll: 



3HriM 



TP 



(46) 



where the proportionality constant has the form A = 
Cr/ P ~ 2 , and C is a dimensionless constant. 

On physical grounds we should expect that p < 3. In 
terms of the correlation length this has the form 



3 dL 
dt 



3HL 



A 



L 2 TP 



(47) 



This may seem very different from the standard term 
given by Eqn. ([7]) , but the difference is actually less than 
it appears. For example, there should be a velocity de- 
pendence in this term (for the obvious reason that if all 
monopoles have zero velocities there will be no annihila- 
tions). This is explicit in the standard term, but implicit 
in the Preskill term. An easy way to see it is to assume 
that monopoles have thermal velocities. In this case one 
can use Eqns. (|23H24|) to put in a phenomenological but 
explicit velocity dependence. We then find 



L 2 TP 



Cv 



(I) 



p-5/2 



(48) 



though we caution that the above assumption of thermal 
velocities is, at best, a crude approximation. In any case, 
we will use the Preskill loss term, Eqn. (|4"T]) . for the lo- 
cal monopoles, and comment of the differences we would 
obtain otherwise. 

As for the specific values of p, Preskill argues for a 
short, transient high-temperature regime where p = 2, 
and a longer low-temperature regime with p = 9/10. This 
is important because if p < 1 annihilations are expected 
to 'turn off' (that is, become unimportant), in which case 
we expect n oc T 3 , which corresponds to L oc t 1 / 2 . On 
the other hand, if p > 1 then annihilations are always 
relevant, and in that case one expects n oc T p+2 , which 
corresponds to I oc £(p+ 2 )/ 6 . (Notice that this analysis 
is for the radiation dominated epoch.) It is therefore 
important to see if we can recover these results using 
the model in the local case. For the global case, a good 
benchmark will be the simulations discussed above. 



V. EVOLUTION OF LOCAL MONOPOLES 

We can now use our newly derived analytic model 
equations to look for scaling solutions for the charac- 
teristic lengthscale and RMS velocity of the monopoles. 
Starting with the case of local monopoles, the evolution 
equations have the general form 



dL 
~dt 



3HL 



L Crf- 



L 2 TP 



and 



dv 
~dt 



(1-v 2 ) 



k L 
7]L 2 d H 



(49) 



(50) 



We can start by finding solutions in Minkowski space- 
time, by setting H = 0. In this case the asymptotic 
scaling solution has the form 



L A cx 



P-2 



-t 



(51) 



and the monopoles will freeze, with the precise scaling 
law for the velocity depending on the behavior of the 
friction length scale. Assuming that we have If = const 
we find 



v oct 



-4/3 



(52) 



while for the arguably more realistic If oc L the freezing 
happens more slowly, 



v oct 



(53) 



In the (unrealistic) frictionless limit If — ► oo the cor- 
relation length still has the same scaling, but velocities 
asymptote to the speed of light, v — > 1. 

For the expanding case, allowing for a generic scale 
factor of the form a oc t x , there are two possible scaling 
laws, which depend both on the value of p and on A. For 
the case 



p< 3 



1 

A 



we will have 



L oc t oc a; 



(54) 



(55) 



this corresponds to the case where energy losses due to 
annihilation are unimportant and the monopoles are sim- 
ply conformally stretched. Note that in the radiation era 
we do recover the result L oc t 1 / 2 , and in that case the 
threshold for this regime is indeed p < 1. This there- 
fore recovers and generalizes the Preskill results. In the 
opposite case 

P > 3 - j (56) 

then annihilations are dynamically important, and the 
scaling law for the correlation length is then 

L oc £( a p +1 )/ 3 ; (57) 
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again we recover the expected Preskill result L oc £(p+ 2 )/ 6 
for the particular case of the radiation era. Note that for 
p > 3-1/A we have (Xp + l)/3 > A: as expected, in this 
regime the evolution is faster than the one above, which 
corresponds to conformal stretching. This difference il- 
lustrates the effect of the annihilations. 

It is also worth pointing out that linear scaling (L cx t) 
will occur for the case Xp — 2, while for Xp > 2 the cor- 
relation length will grow supcrluminally. This last case 
corresponds (at a phenomenological level) to the situa- 
tion where annihilations are so efficient that on average 
there will eventually be less than one monopolc per Hub- 
ble volume, so the monopoles effectively disappear. Note 
that since we physically expect p < 3, then Xp = 2 corre- 
sponds to A > 2/3. Hence linear scaling can occur in the 
matter era but not in the radiation era. On the other 
hand, super luminal scaling requires A > 2/3 and so it 
can't occur in cither epoch — this is a simple manifesta- 
tion of the monopole problem in standard cosmology. 

Interestingly, in both regimes the scaling law for the 
velocities is the same, namely 



!)cxr A «fl _1 cxT. 



(58) 



This is a nice and simple result, and it disproves the 
naive expectation that the monopoles should move with 
thermal velocities (which would correspond to v cx VT). 

The above solutions hold for decelerating universes 
(with < A < 1) but also for power-law inflating uni- 
verses (that is, with A > 1). On the other hand, in de 
Sitter space (with a cx e Ht ) we have 



or 



L cx a, p < 3 



Loca p/3 , p>3. 



(59) 



(60) 



although we expect that the latter behavior for p is phys- 
ically unrealistic. For the velocity we still have 



v cx a 



(61) 



Thus in an inflating universe the monopoles freeze and 
are conformally stretched, being pushed outside the hori- 
zon. After the end of inflation there will be much less 
than one monopole per horizon and their velocities will 
be infinitesimal, so they will keep being conformally 
stretched until they re-enter the horizon. So in order to 
solve the monopole problem one needs sufficient e-folds 
of inflation to ensure that the monopoles have not yet 
rc-cntcrcd. 



VI. EVOLUTION OF GLOBAL MONOPOLES 

A similar analysis can now be done for the global case. 
We shall see that the different force and energy loss terms 
will lead to very different scaling laws. This case is also 
interesting because previous numerical simulations exist 



against which we can compare our results. Although the 
simulations have been done several years ago and have 
very low resolution by today's standards, we shall see 
that the results of the comparison are very encouraging. 



A. Analytic Model 

In this case the evolution equations will have the gen- 
eral form 



3— = 3HL + v - — V cv 

dt t d 



and 



k f±\ 3/2 _v_ 
L \d H J i d 



(62) 



(63) 



Again we can start with the Minkowski space-time 
case. In the unrealistic frictionless limit if — > oo we 
now have asymptotically 



L=-ct 



v = 1 



(64) 



(65) 



so global monopoles will become ultra-relativistic. The 
case of a constant friction length scale is not relevant for 
global monopoles. Due to the linear divergence of their 
masses, a more realistic situation in Minkowski space 
time would be that of the friction length scale being pro- 
portional to the correlation length itself, If cx L. In that 
case we find the following scaling law 



L = et = -v a (v + c)t 



fee 3 / 2 



const. 



(66) 



(67) 



This behavior is to be contrasted with the case of local 
monopoles, whose velocities always approach zero (ex- 
cept in the unrealistic frictionless case). Note that in 
principle any value of the velocity is a possible solution, 
including the limit v = 1. An interesting question is 
whether the friction will make the monopole velocities 
stabilize at some fixed value (and if so, how small this is) 
or if they will still become arbitrarily close to the speed 
of light. Incidentally, notice that assuming if cx t a in 
Minkowski space, requiring a linear scaling L cx t im- 
plies v = const and o = 1. In other words, no other 
non-trivial behavior of the friction length would lead to 
linear scaling for the network. 

Now let us consider the general expanding case, again 
for a generic expansion law a oc t A . Here, just as in the 
standard case of local strings, the only possible scaling 
law is linear scaling 



L = et 



(68) 
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v = vq = const., (69) 

which at least qualitatively is in agreement with the ex- 
isting numerical simulations. Interestingly, just as in 
the Minkowski case there are two branches of the solu- 
tion: the velocities may or may not be ultra-relativistic! 
Firstly, the ultra-relativistic scaling regime will have the 
following scaling parameters 

vq = 1 (70) 



other scaling solutions exist. The only change is that in 
the scaling coefficients (e and Vq above), we would need 
to interpret the parameter A as having a renormalized 
value, instead of the value given by the expansion rate. 

We can also find scaling solutions in inflating universes. 
Here there is a unique solution, both for power-law infla- 
tion and for the de Sitter case, namely 

Laa (77) 



3-4A' 



(71) 



note that this can only hold for A < 3/4, but is in princi- 
ple allowed both in the radiation and in the matter eras. 
Secondly, the more standard (sub-luminal) scaling regime 
will be characterized by the following scaling parameters 



cv 



3(1 - A) - \v 2 



Xv = fc(l - A) 



3/2 1/2. 



(72) 



(73) 



These relations could be solved explicitly for e and Vq, 
but the corresponding expressions would not be too illu- 
minating. However, simplified and physically suggestive 
solutions can be displayed for both limits of the expan- 
sion power A. In the limit A — > 0, we have 



1 



-CV . 



(74) 



Not surprisingly, this is similar to the Minkowski space- 
time scaling we discussed above. On the other hand, in 
the limit A — > 1, we find 



v = r i ck 2 (l-\) 2 



(75) 



Here the scaling velocity becomes arbitrarily small (v — > 
0) and L oc a cx t so asymptotically this is a conformal 
stretching regime. 

Unlike the ultra-relativistic branch, this non-luminal 
branch can exist for any A (that is, for any expansion 
law), though note that there is a constraint on the scaling 
value of the velocity 



1 



< < 3( j - 1) 



(76) 



this is trivial for A < 3/4 (meaning that in such cases any 
scaling value for the velocity is allowed in principle), but 
restrictive for faster expansion rates. Note in particular 
that it agrees with the above finding that v = 1 is only 
allowed for A < 3/4. 

In passing we also note that the linear scaling solution 
L oc t, v = const will also hold if we consider the case 
where the friction length varies as if oc L instead of the 
usual scaling with temperature. Even in that case no 



woe a -1 ; (78) 

hence the solution is the same as in the local case, and 
the number of e-folds of inflation required to keep the 
monopoles outside the horizon by the present day should 
be the same in both cases. 



B. Comparing to Simulations 

Our results suggest that it would be very interesting to 
carry out high-resolution numerical simulations of global 
monopoles with a range of different expansion rates (say, 
radiation, matter and a value of A > 3/4) in order to 
check these solutions and provide a good calibration for 
the model. In the meantime, however, we can use the 
results of Bennett and Rhic [18] and of Yamaguchi Oil 



for the correlation length, and also the latter's for the 
velocities, in order to make some simple comparisons. 

We start by translating these results into our scaling 
parameter e, finding 

e r ~ 1.32, e m ~ 1.89 (79) 

for the simulations of Bennett and Rhie, and 



1.32, 



1.59 



(80) 



for Yamaguchi's. Notice the remarkable agreement in 
the radiation era; even in the matter era the difference 
is small considering the relatively low resolution and dy- 
namic range of the simulations. Now, since Yamaguchi's 
measurement is consistent with luminal velocities, let us 
assume that we are in the v = 1 branch, and solve for 
the energy loss parameter c. We then find 



c r ~ 1.32, 



0.63 



for the simulations of Bennett & Rhie, and 



1.32, 



0.53; 



(81) 



(82) 



notice that there is a factor of two difference between 
the values of the parameter in the radiation and matter 
eras, while we would expect to find similar values in both 
epochs if the model is broadly correct and the parameter 
c is a constant (or nearly so). On the other hand, if we 
assume that we are in the sub-luminal branch (using in 
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both cases Yamaguchi's values for the velocities, since 
Bennett and Rhie doesn't provide a measurement), 



Cf 1.32, Cm 1-35 
for the simulations of Bennett and Rhie, and 



1.32, 



1.15; 



(83) 



(84) 



here we can claim a very good agreement, considering the 
resolution of the simulations in question and the error 
bars that each of them has. Also in this branch we can 
compare the scaling values of the velocities; according to 
the model we expect the ratio between the matter and 
radiation era scaling values to be 



(85) 



while Yamaguchi finds v m /v r ~ 0.8; again given the error 
bars we would argue that the agreement is encouraging. 



VII. REVISITING OUR ASSUMPTIONS 

As we discussed in Sect. IIV1 most of the subtlety in 
the analytic description of these defect networks rests 
in the way the energy losses, the inter-monopole forces 
and their suppression factors are incorporated into the 
model. In this section we will revisit these assumptions 
and discuss the extent to which they influence the results 
we presented so far. 

Let us start by commenting on the role of the anni- 
hilation term, by considering what would happen if we 
had used PreskuTs annihilation term — Eqn. (|4"Tj) , which 
holds for local monopoles — for the global monopole anal- 
ysis of Sect. I VII For scaling laws of the form L a 
d a t^, we find the following scaling powers for the cor- 
relation length 



= A, p<3- 



and 



Xp+ 1 



p>3-~ 



(86) 



(87) 



(again annihilations are unimportant in the former case 
but not in the latter). Not surprisingly the scaling laws 
for L would be similar to those we found in the local case 
(Sect. [V}. On the other hand, for the velocities we have 



and 



/3 = -\, A<l/3 



0=-S-, A>l/3. 



(89) 



Again the first of the above existed in the local case (for 
all values of A) , while the second branch is new. So there 



would in general be no linear scaling solution; indeed the 
only possibility of having a = 1 is for the special case 
Xp = 2 (with p > 2). But in any case the velocities 
will become arbitrarily small, so as expected this form 
of annihilation term seems clearly ruled out even by the 
existing, low-resolution numerical simulations. 

A more important issue is the suppression factor on the 
forces driving the monopoles, to account for the partial 
cancellation due to the presence of many monopoles and 
antimonopolcs in a given Hubble volume. In Sect. IIVI 
we argued that this could be described as a 1/y/N effect, 
but with the number N of defects in a Hubble volume 
calculated differently for the local and global cases. Here 
we will generalize this assumption, by assuming that N 
is generically given by 



Y s 



L 



(90) 



recall that in the above analysis it was assumed that s = 
for the global case but s = 1 for the local case. 

Our analysis in Sects. FVl and |VT1 can now be repeated 
for this generic suppression term. The general outcome is 
that our previous analysis is fairly robust. Most scaling 
laws are not affected or receive s-dependent corrections 
only in the pre- factors. Corrections in the scaling expo- 
nents appear only for the velocities, not the correlation 
lengths. On the other hand, the emergence of genuinely 
new scaling regimes is only possible for (arguably) unre- 
alistic values of the parameter s. We now discuss these 
changes separately for the local and global cases. 



A. Local case 

The evolution equations now take the form 

L Crf- 2 



dL 

3— = 3HL- 
dt 



id 



L 2 TP 



and 



dv 
11 



(I-- 2 ) 



— (-) 



(3- S )/2 



V 

~l 



(91) 



(92) 



Recall that our expectation, based on the results of Ein- 
horn et at. [38[ is that s = 1: the suppression factor 
describing the partial cancellation of the force term is 
weaker due to monopole-antimonopole anticorrelations, 
which can be reinterpreted as an reduction in the effec- 
tive number of neighbours. 

Starting with Minkowski case (H = 0), there are no 
changes in the asymptotic case If — > oo. For if = const 
the scaling law for the lengthscale L (given in Eqn. (fBTj) ') 
also remains unchanged, for any values s < 7/2, On the 
other hand, the power law of the scaling of the velocities 
(v oc t") does have an s-dependence 



(93) 
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for s = 1 we therefore recover our previous result in Eqn. 
(f52"| . (3 = —4/3. Something analogous happens for the 
arguably more realistic case if cx L: the scaling law for 
L is unchanged for s < 3, while the power law of the 
velocities becomes 

P = L ^\ (94) 

again for s — 1 we recover Eqn. (f53"| . /3 = — 1. In this 
case there is in principle a new scaling branch for s > 3 
(even though we emphasize that such a suppression factor 
is physically unrealistic). In this case we would have 
(3 = — 1/s and the characteristic lengthscale would also 
scale differently: L cx t a with a = 1 — (2/s). 

For the case of the expanding universe, we again find 
that the scaling law for the lengthscale remains un- 
changed, being given by Eqns. ([55| or (|57|) as before, 
depending on the value of the Preskill factor p. For the 
velocities there are two possibilities. If s < 1 the scal- 
ing law is again unchanged, being given by (3 = — A (for 
a cx t x ) as in Eqn. On the other hand, for s > 1 the 

velocities will scale differently, with the scaling exponent 
being given by 

/3 = -f(s + l) + ~( S -l), (95) 

where again a is the scaling exponent for L given by 
Eqns. ([55]) or (|57|) . Careful measurements of the scaling 
laws for L and v in numerical simulations with several 
expansion rates can therefore be used to not only test the 
analytic model but also to obtain indirect information on 
the inter-monopole forces. 



which trivially reduces to Eqn. (|67)) for s — 0. (In partic- 
ular, monopole velocities arbitrarily close to the speed of 
light are still possible in principle.) Interestingly, for the 
case s = 1 only, there is a second scaling solution: writing 
L cx t a and v cx t 13 as usual, this solution is characterized 
by 

a = 1 + (3 = -ck , (99) 

subject to the constraints that a < 1 and (3 < which 
trivially imply ck < 3. Unlike the generic linear scaling 
regime, here the characteristic lengthscale grows more 
slowly and the velocities become arbitrarily small instead 
of being constant. 

Something similar occurs in the case of the expanding 
universe. The linear scaling solution still exists for all 
values of s (including the branch with v = 1), the only 
difference being an s-dependence in the relation between 
the pre-factors for the L and v scaling laws, formerly 
given by Eqn. (|73p and now becoming 

Xv Q = k(l - A^-sVV 1 ^)/ 2 , (100) 

which trivially reduces to the previous result when s = 0. 

In this case a second, non-relativistic scaling solution 
appears only for the case s = — 1. The physical viability 
(if any) and interpretation of such a term is somewhat 
unclear: this would imply that the suppression factor is 
stronger than the naive l/y/~N effect, or in other words 
each monopole has a larger effective number of neigh- 
bours. Just as in the Minkowski case, the new scaling 
solution is characterized by 



B. Global case 



a = l + (3. 



(101) 



In this case the evolution equations will have the 
generic form 



dL o L 

3— = 3HL + v — — h cv 

at £ d 



and 



dv 
~dt 



fc/iX (3- S )/2 v 



L\d 



(96) 



(97) 



where thus far we have assumed s = 0, corresponding to 
the expectations that in the absence of any correlation or 
anticorrelations the suppression factor should be a simple 
1/y/N effect. 

Again starting with the Minkowski space case (H = 0), 
the case if — > oo is still unchanged. For the case if cx L 
the linear scaling solution L = et, v — const given by 
Eqns. (|6"6"|) and ([6T)l still exists for all values of s, the 
only difference being that the relation between the pre- 
factors becomes s-dependent 



v = fce (3-s)/2 , 



(98) 



subject to a < 1 and (3 < 0, but the specific values of the 
scaling coefficients now depend on the expansion rate, 
being given by 



1 



03 + A)(l + /3-A) - -(1- X) 2 ck 



(102) 



and subject to the further constraints that both terms in 
brackets on the left-hand side of the above equation are 
positive. 

Although this scaling solution is clearly related to the 
one we discussed previously in Minkowski space, one 
can't recover the Minkowski one simply by taking the 
limit A — > 0: this is due to the fact that the behavior of 
the damping terms is different in the two cases (It oc L 
for the Minkowski case versus Hubble damping for the 
expanding case). This is also the reason why different 
values of s are needed in the two cases for the solution 
to exist. The two new scaling solutions have the distin- 
guishing feature of velocities decreasing as the network 
evolves and becoming arbitrarily small (as opposed to be- 
ing constant) , which again highlights the need for careful 
measurements of monopole velocities. 
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VIII. SUMMARY AND OUTLOOK 

We have proposed a velocity-dependent one-scale phe- 
nomenological description of the evolution of a network 
of pointlike defects -monopoles-, both in the global and 
local (magnetic) cases. In both cases we recover scal- 
ing solutions previously found in the literature and make 
predictions for new scaling regimes that have not been 
considered before. 

In the global monopole case there are three main 
assumptions. First, since the force between global 
monopoles is approximately independent of distance, a 
monopole interacts with every other monopole and anti- 
monopole within its Hubble volume and this reduces the 
effective force it feels. Second, since its mass grows lin- 
early with distance the acceleration has the same 1/L(t) 
dependence as for local cosmic strings. Finally, since 
radiation losses occur mainly within bound monopolc- 
antimonopole pairs, and are therefore decoupled from the 
rest of the network in a first approximation, they can be 
modelled by a single energy-loss parameter decribing the 
efficiency of pair formation. Comparison with existing 
numerical simulations [H, [l!| supports this overall pic- 
ture and we recover scaling solutions with v < 1 previ- 
ously reported in the literature. The rms velocity of the 
monopoles is not determined by the model. On the other 
hand, the expected v — 1 scaling solution does not agree 
so well with the existing simulations. However the calcu- 
lation of monopole velocities from numerical simulations 
has large uncertainties so it would be important to revisit 
this point numerically. 

In the case of magnetic monopoles we build on 
Preskill's results [25[ for the annihilation rate. The sup- 
pression of the Coulomb force between monopoles due to 
the presence of many neighbours is further corrected to 
account for anticorrelations (a monopole's neighbour is 
more likely to be an antimonopole). 

Mindful of the fact that our modeling relies on a few 
key assumptions, we have also discussed how our results 
change if some of these assumptions are relaxed. We have 
found that our scaling solutions are usually quite robust, 
though more so for the defect correlation lengths than for 
their velocities. In particular there can be changes if the 
suppression factors for the forces driving the monopoles 
are different. A context where this may possibly happen 
is that of condensed matter systems such as Helium or 
liquid crystals. The behavior of monopole-like defects in 
these systems deserves further study. 



It would be very interesting to extend this analysis 
to the case of hybrid semilocal networks (of monopoles 
connected by strings), for which ref. [40] found some 
evidence of scaling. Numerical simulation is computa- 
tionally very demanding for these systems, because one 
cannot rely on any thin string or sigma model approxi- 
mations. 

In semilocal networks, the long strings and loops be- 
have like local cosmic strings but string segments are also 
possible and they can dominate the dynamics (see [4l| for 
a review). At formation the network consists only of seg- 
ments [42J, with an exponential distribution in length. 
But the ends of segments have long-range interactions 
that make them behave much like the global monopoles 
studied here. Once formed, they will interact and anni- 
hilate with other segment ends even at long distances, 
causing some segments to grow into long strings while 
others will collapse and disappear. 

This long-range interaction provides a crucial differ- 
ence with the better studied case of hybrid networks of 
monopoles connected by strings in which the strings con- 
fine the magnetic flux of the monopoles [2]. In these net- 
works the ends of the strings are very light and the dy- 
namics is dominated by the tendency of the segments to 
collapse. In the semilocal case, which of these two effects 
dominates the network evolution is controlled by the ra- 
tio of the scalar and gauge couplings (see also [43| for a 
study of the electroweak case). For large scalar quartic 
coupling the strings disappear, leaving behind a scaling 
network of texture-like structures. But if the coupling 
is small the simulations hint at a scaling network made 
of local strings plus a few open segments per horizon 
volume. If such structures are indeed formed after e.g. 
brane inflation, an extension of the present study would 
be the right tool to confirm or rule out scaling behavior. 
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